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I> ! Abstract 

00 . 

In this talk CP violation in the supersymmetric models, and especially in B- 
■^j- ' decays is discussed. We review our analysis of the supersymmetric contributions to 

the mixing CP asymmetries of B — > <pK$ and B — > rj ' K$ processes. Both gluino 
and chargino exchanges are considered in a model independent way by using the 
(— 1 , mass insertion approximation method. The QCD factorization method is used, and 

Qlh! parametrization of this method in terms of Wilson coefficients is presented in both 

Q_i! decay modes. Correlations between the CP asymmetries of these processes and the 

direct CP asymmetry in b — > sj decay are shown. 

> : 
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1 Introduction 



In the Standard Model (SM), the CP violation is due to the misalignment of the mass 
and charged current interaction eigenstates. This misalignment is represented in the CKM 
matrix Vqkm [1], present in the charged current interaction Lagrangian, 
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In the Wolfenstein parametrization Vckm is given by 
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where the Cabibbo mixing angle A = 0.22. The CKM matrix is unitary, Vq KM Vckm = 
1 = VckmVq KM . The unitarity conditions provide strong constraints for CP violation 
in the Standard Model. It is, however, well known that the amount of CP violation in 
the Standard Model is not enough to account for the generation of the matter-antimatter 
asymmetry in the universe. Thus new sources for CP violation are expected from beyond 
the Standard Model scenarios. E.g. in supersymmetric models a large number of new 
phases emerge. These phases are strongly constrained by electric dipole moments. 

The unitarity constraints can be represented as unitarity triangles, for which the length 
of the sides correspond to the products of elements in the CKM matrix due to the unitarity 
conditions. In the Standard Model, the angle (3 in the unitary triangle [2] , can be measured 
from B meson decays. The golden mode B° — > J/ ipKs is dominated by tree contribution 
and measurement of the CP asymmetries very accurately gives the (5 angle. 

The dominant part of the decay amplitudes for B° — > <j)Ks, rj'Ks is assumed to come 
from the gluonic penguin, but some contribution from the tree level b — > uus decay is 
expected. The \<f>) is almost pure \ss) and consequently this decay mode corresponds 
also accurately, up to terms of orders 0(X 2 ), to sin 2/3 in the SM [3]. The b — > uus tree 
level contribution to B d — > rj'K was estimated in [4]. It was found that the tree level 
amplitude is less than 2% of the gluonic penguin amplitude. Thus also in this mode one 
measures the angle (3 with a good precision in the SM. Therefore, it is expected that NP 
contributions to the CP asymmetries in B° — > <pKs, rj'Ks decays are more significant 
than in B° — > J/ipK s and can compete with the SM one. 

B physics is a natural framework to test the beyond the Standard Model CP violation 
effects. It is clear that ultimately one needs to test the Standard Model with three 
generations and that large statistics is needed to achieve the small effects of CP violation. 
Flow of interesting data has been provided in recent years by the B-factories. 

In this talk we will not assume any particular model for CP violation, but consider 
a general SUSY model. The talk is based on papers [5,6], where a comparative study 
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of SUSY contributions from chargino and gluino to B — > <pK and B — > rj'K processes 
in naive factorization (NF) and QCD factorization (QCDF) approaches is done. We also 
analyzed in [6] the branching ratios of these decays and investigated their correlations 
with CP asymmetries. The correlations between CP asymmetries of these processes and 
the direct CP asymmetry in b — > S7 decay [7] is also discussed in [6]. 

In the analysis the mass insertion method (MIA) [8] is used. Denoting by A q AB the 
off-diagonal terms in the sfermion (q = u, d) mass matrices for the up and down, respec- 
tively, where A, B indicate chirality couplings to fermions A, B — (L, R), the A-B squark 
propagator can be expanded as 

«j£> = i (kH - fn'l - A^)" 1 K + + O(A'), (3) 

where q = u,d selects up or down sector, respectively, a,b — (1,2,3) are flavor indices, 
1 is the unit matrix, and rh is the average squark mass. As we will see in the following, 
it is convenient to parametrize this expansion in terms of the dimensionless quantity 
(5%) ab ee (A%) ab /m*. 

New physics (NP) could in principle affect the B meson decay by means of a new 
source of CP violating phase in the corresponding amplitude. In general this phase is 
different from the corresponding SM one. If so, then deviations on CP asymmetries from 
SM expectations can be sizeable, depending on the relative magnitude of SM and NP 
amplitudes. For instance, in the SM the B — > <j)K s decay amplitude is generated at 
one loop and therefore it is very sensitive to NP contributions. In this respect, SUSY 
models with non minimal flavor structure and new CP violating phases in the squark 
mass matrices, can easily generate large deviations in the B — > <pKs asymmetry. 

The time dependent CP asymmetry for B — > <pKs can be described by 

T(B°(t) -> <f)K s ) - T(B(t) -> <j)K s ) „ . n , 

a<pKs (t) = -A^J <LM ^l = C^ K8 co6AM Bd t + S 4Kg8 mAM Bd t, 

T(B (t) ->■ (pK s ) + T(B(t) -> (pK s ) 

(4) 

where C^k s an d S<\>k s represent the direct and the mixing CP asymmetry, respectively 
and they are given by 



\p{<t>K s )\ 2 -l 2/m 



(5) 
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The parameter p((f>Ks) is defined by 

p{<pKs) = A^KiY (6) 

where A((f>Ks) and A(<fiKs) are the decay amplitudes of B° and B° mesons, respectively. 
Here, the mixing parameters p and q are defined by \Bi) = p\B) + q\B ), \B 2 ) = 
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p\B) — q\B ) where |i?i(2)) are mass eigenstates of B meson. The ratio of the mixing 
parameters is given by 

1 - r- 2id * Vt * bVtd (7) 

n~ V,V*' ^ ' 

V v tb v td 

where 9 d represent any SUSY contribution to the B — B° mixing angle. Finally, the above 
amplitudes can be written in terms of the matrix element of the AB = 1 transition as 

A^Ks) = (tPKslH^lB*), A{4>K S ) = (<t>K s \ (H^f \B°). (8) 

Results by BaBar and Belle collaboration have been announced in [9,10]. The ex- 
perimental value of the indirect CP asymmetry parameter for B° — > J/ipKs is given by 
[9, 10] 

Sj/^k s = 0.726 ± 0.037, (9) 

which agrees quite well with the SM prediction 0.715ig;o45 [11]. Results on the corre- 
sponding sin 2/3 extracted for B° — > 4>K S process is as follows [9, 10] 

S 4>Ks = 0.50±0.25t°;°J (BaBar), 

= 0.06 ±0.33 ±0.09 (Belle), (10) 

where the first errors are statistical and the second systematic. As we can see from Eq.(10), 
the relative central values are different. BaBar results [9] are more compatible with SM 
predictions, while Belle measurements [10] show a deviation from the cc measurements of 
about 2a. Moreover, the average S^Ks = 0.34 ±0.20 displays 1.7cr deviation from Eq.(9). 

Furthermore, the most recent measured CP asymmetry in the B° — > r/'Ks decay is 
found by BaBar [9] and Belle [10] collaborations as 

S^Ks = 0.27 ±0.14 ±0.03 (BaBar) 

= 0.65 ±0.18 ±0.04 (Belle), (11) 

with an average S v 'k s = 0-41 ± 0.11, which shows a 2.5<r discrepancy from Eq. (9). 

It is interesting to note that the results on s-penguin modes from both experiments 
differ from the value extracted from the cc mode (J/ip), BaBar by 2.7cr and Belle by 2.4<j 
[9,10]. At the same time the experiments agree with each other, and even the central 
values are quite close: 

0.42 ±0.10 BaBar, 0.43^ Belle. 

On the other hand, the experimental measurements of the branching ratios of B° — > 
<pK° and B° -> r]'K° at BaBar [12], Belle [13], and CLEO [14] lead to the following 
averaged results [15] : 

BR(B° -> (j)K°) = (8.3±i;g) x 10" 6 , (12) 
BR(B° -> V 'K°) = (65.21™) x 10" 6 . (13) 
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From theoretical side, the SM predictions for BR(B — > <f>K) are in good agreement with 
Eq.(12), while showing a large discrepancy, being experimentally two to five times larger, 
for BR(B — > rj'K) in Eq.(13) [16]. This discrepancy is not new and it has created a 
growing interest in the subject. However, since it is observed only in B — > rj'K process, 
mechanisms based on the peculiar structure of rj' meson, such as intrinsic charm and 
gluonium content, have been investigated to solve the puzzle. Correlations with branching 
ratios have been discussed in [6]. 



2 SUSY contributions to B -> (j>(r/)K decay 

We first consider the supersymmetric effect in the non-leptonic AB = 1 processes. Such 
an effect could be a probe for any testable SUSY implications in CP violating experiments. 
The most general effective Hamiltonian H^ B=1 for the non-leptonic AB = 1 processes 
can be expressed via the Operator Product Expansion (OPE) as [17] 



ttAB=1 
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where \ p = V pb V* s , with V pb the unitary CKM matrix elements satisfying the unitarity 
triangle relation A t + A M + A c = 0, and Cj = Cj(/ib) are the Wilson coefficients at low 
energy scale fj, b ~ 0(m b ). The basis Q, L = Qi(fib) is given by the relevant local operators 
renormalized at the same scale /j b , namely 
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Here a and /3 stand for color indices, and t^g for the SU(3) C color matrices, a 



(15) 
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|i[7^, 7"]. Moreover, e g are quark electric charges in unity of e, (??)y±A = 97/i(l ± 75)9, 
and q runs over -u, d, s, and 6 quark labels. In the SM only the first part of right hand side 
of Eq.(14) (inside first curly brackets) containing operators Qi will contribute, where Q p l 2 
refer to the current-current operators, Q3-6 to the QCD penguin operators, and Q7-10 to 
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the electroweak penguin operators, while Q-j^ and Q 8g are the magnetic and the chromo- 
magnetic dipole operators, respectively. In addition, operators Qi = (/•*&) are obtained 
from Qi by the chirality exchange (qiq2)v±A — > (qi^vtA- Notice that in the SM the 
coefficients Ci identically vanish due to the V-A structure of charged weak currents, while 
in MSSM they can receive contributions from both chargino and gluino exchanges [18, 19]. 

As mentioned, we calculated in [5] the chargino contribution to the Wilson coefficients 
in the MIA approximation. In MIA framework, one chooses a basis (called super-CKM 
basis) where the couplings of fermions and sfermions to neutral gaugino fields are flavor 
diagonal. In this basis, the interacting Lagrangian involving charginos is given by 

Asx+ = -9 EE (Vki K* ba d a L (xtY u\ - U* k2 (Y d di ^.K + ) ab 4 (xtr u\ 

k a,b 

-V k2 (K.Y u d ^) ab dl(xtru b R ), (16) 

where q^i — |(1 ± 75)?) and contraction of color and Dirac indices is understood. Here 
Y^™ g are the diagonal Yukawa matrices, and K stands for the CKM matrix. The indices 
a, b and k label flavor and chargino mass eigenstates, respectively, and V, U are the 
chargino mixing matrices defined by 

WM^V-' = diagK + ,m- + ), and M# I - ^ Q v^^sin/? \ 

where M 2 is the weak gaugino mass, \x is the supersymmetric Higgs mixing term, and tan f3 
is the ratio of the vacuum expectation value (VEV) of the up-type Higgs to the VEV of the 
down-type Higgs 1 . As one can see from Eq.(16), the higgsino couplings are suppressed by 
Yukawas of the light quarks, and therefore they are negligible, except for the stop-bottom 
interaction which is directly enhanced by the top Yukawa (Y t ). In our analysis we neglect 
the higgsino contributions proportional to the Yukawa couplings of light quarks with the 
exception of the bottom Yukawa Y b , since its effect could be enhanced by large tan/5. 
However, it is easy to show that this vertex cannot affect dimension six operators of the 
effective Hamiltonian for AB = 1 transitions (operators <5i=i-io in Eq.(14)) and only 
interactions involving left down quarks will contribute. On the contrary, contributions 
proportional to bottom Yukawa Y b enter in the Wilson coefficients of dipole operators 
(C 77 , C 8g ) due to the chirality flip of b — > S7 and b — > sg transitions. 

The calculation of B — > 4>(r]')K decays involves the evaluation of the hadronic matrix 
elements of related operators in the effective Hamiltonian, which is the most uncertain 
part of this calculation. In the limit in which m b ^> Aq C d and neglecting QCD corrections 
in a s , the hadronic matrix elements of B meson decays in two mesons can be factorized, 
for example for B — > M1M2, in the form 

(M.M^Q.IB ) = (M^B ) x (M 2 \j 2 \0) (18) 

1 This tan/3 should not be confused with the angle /3 of the unitarity triangle. 
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where M 12 indicates two generic mesons, Qi are local four fermion operators of the effec- 
tive Hamiltonian in Eq.(14), and represent bilinear quark currents. Then, the final 
results can be usually parametrized by the product of the decay constants and the transi- 
tion form factors. This approach is known as naive factorization (NF) [20,21]. In QCDF 
the hadronic matrix element for B — > MK with M = (f>,rf in the heavy quark limit 
m b ^> Aqcd can be written as [22] 

(MK\Qi\B) qcdf — {MK\Qi\B) NF- 

where (MK\Qi\B)^p denotes the NF results. The second and third term in the bracket 
represent the radiative corrections in 0:5 and AQCD/ m b expansions, respectively. Notice 
that, even though at higher order in a s the simple factorization is broken, these corrections 
can be calculated systematically in terms of short- distance coefficients and meson light- 
cone distribution functions. 

In the QCD factorization method [22,23], the decay amplitudes of B — > 4>(rj')K can 
be expressed as 

A (B -> <f>(r]')K) = A f (B^ (f>(r]')K) + A a (B -> 4>{r]')K) , (20) 

where 

10 

A* (B - <ftTf)K) = ^EE V Pb V; s at^immmNF, (21) 

V " p=u,c 1=1 

and 

(~i 10 
A a (B - <P( V ')K) = ^/b/kU E E W bt W) . (22) 

p=u,ci=l 

The first term A* (B — > (f>(j]')K) includes vertex corrections, penguin corrections and hard 
spectator scattering contributions which are involved in the parameters af^ v \ The other 
term A a (B — > (f)(r]')K) includes the weak annihilation contributions which are absorbed 
in the parameters bp v \ However, these contributions contain infrared divergences, and 
the subtractions of these divergences are usually parametrized as [23] 

[ l * ^X B * S J (23) 

Jo x v ) \AqcdJ 

where p H ,A are free parameters expected to be of order of p H ,A — 0{\), and (/>h,a £ [0, 2n]. 
As already discussed in Ref.[23], if one does not require fine tuning of the annihilation 
phase (f>A, the pa parameter gets an upper bound from measurements on branching ratios, 
which is of order of pa ^ 2. Clearly, large values of pa are still possible, but in this 
case strong fine tuning in the phase (f>A is required. However, assumptions of very large 
values of Ph,a, which implicitly means large contributions from hard scattering and weak 
annihilation diagrams, seem to be quite unrealistic. In [6] we assumed p < 2. 
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Figure 1: S$k s as a function of arg[( 5^) ] (left) and argH^^j ] (right) with gluino con- 

\ / 23 ^ * 23 

tribution of one mass insertion. The region inside the two horizontal lines corresponds to the 
allowed experimental region at la level. 

3 CP asymmetry in B —> <j>Ks and in B —> rj'Ks 

In order to simplify our analysis, it is useful to parametrize the SUSY effects by introducing 
the ratio of SM and SUSY amplitudes as follows 



4>Ks 



and analogously for the rj'Ks decay mode 

/ j\ SVSY \ 
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where Ri stands for the corresponding absolute values of | a ^m | , the angles 6s, ^ are 
the corresponding SUSY CP violating phase, and 8s, v > = 8^, — 8^5* parametrize the 
strong (CP conserving) phases. In this case, the mixing CP asymmetry Ssk s in Eq.(4) 
takes the following form 

sin 2(3 + 2R 6 cos 5s sin(0^ + 2(3) + Rl sin(2^ + 2(3) 



S, 



<t>K s 



1 + 2R 6 cos <L cos 6s + Ra 



(27) 



and analogously for B — > rj'Ks 
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sin 2(3 + 2R V , cos 8^ sin(^/ + 2(3) + R\, sin(20^ + 2(3) 
1 + 2R T] < cos 8 V > cos 9 v i + R 2 , 



(28) 
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Arg[(5 a LL ) 23 ] 




Arg[(5 LR ) 23 ] 



Figure 2: As in Fig. 1, but for S^Ks- 



Our numerical results for the gluino contributions to CP asymmetry S$k s are P re ~ 
sented in Fig. 1, and to CP asymmetry S v 'k s are presented in Fig. 2. In all the plots, 
regions inside the horizontal lines indicate the allowed 2a experimental range. In the plots 
only one mass insertion per time is taken active, in particular this means that we scanned 
over I < 1 an d \ {^lr) I < 1- Then, Sq^Ks * s plotted versus 9$, which in the 

\ / 23 ^ * 23 

case of one dominant mass insertion should be identified here as 9$ = a.rg[(5 AB )ij}. 

We have scanned over the relevant SUSY parameter space, in this case the average 
squark mass rh and gluino mass m g , assuming SM central values [24]. Moreover, we 
require that the SUSY spectra satisfy the present experimental lower mass bounds [24]. 
In particular, rrig > 200 GeV, rh > 300 GeV. In addition, we impose that the branching 
ratio (BR) of b — > 57 and the B — B mixing are satisfied at 95% C.L. [25], namely 
2 x 10~ 4 < BR(b -> S7) < 4.5 x 10~ 4 . Then, the allowed ranges for | (df L ) \ and | (<Jf fl ) | 

V / 23 * f 23 

are obtained by taken into account the above constraints on b — > sj and B — B mixing. 
We have also scanned over the full range of the parameters pa,h and 4>a,h in Xa an d Xh, 
respectively, as defined in Eq.(23). 

The chargino effects to S^k s and S v 'k s [5, 6] are summarized in Fig. 3 and Fig. 4, 
where S^)k s * s plotted versus the argument of the relevant chargino mass insertions, 
namely (Sl L ) 32 and (5% L ) 32 . 

As in the gluino dominated scenario, we have scanned over the relevant SUSY param- 
eter space, in particular, the average squark mass rh, weak gaugino mass M 2 , the \i term, 
and the light right stop mass vn~ tR . Also tan/3 = 40 has been assumed and we take into ac- 
count the present experimental bounds on SUSY spectra, in particular rh > 300 GeV, the 
lightest chargino mass M x > 90 GeV, and rhi R > 150 GeV. As in the gluino case, we scan 
over the real and imaginary part of the mass insertions (<^ L ) 32 and {S^ L ) 32 , by considering 
the constraints on BR(6 — > 57 ) and B — B mixing at 95% C.L.. The b — > 57 constraints 
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Figure 3: As in Fig. 1, but for S$k s as a function of a^gKS^)^] (left) and arg[(<^ L ) 32 ] (right) 
with chargino contribution of one mass insertion . 
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Figure 4: As in Fig. 3, but for S v >k s - 



impose stringent bounds on (S^ L ) 32 , specially at large tan/5 [5]. Finally, as in the other 
plots, we scanned over the QCDF free parameters pa,h < 2 and < 4>a,h < 27r. 

The reason why extensive regions of negative values of S^Ks are excluded here, is 
only due to the b — > 57 constraints [5]. Indeed, the inclusion of (5^)32 mass insertion 
can generate large and negative values of S^x s , by means of chargino contributions to 
chromo-magnetic operator Q 8g which are enhanced by terms of order m^jm^. However, 
contrary to the gluino scenario, the ratio |Cs s /CV 7 | is not enhanced by color factors and 
large contributions to Cs g leave unavoidably to the breaking of b — > S7 constraints. 

We plot in Figs. 5 the correlations between S^Ks versus S v >k s for both chargino and 
gluino in QCDF. For illustrative purposes, in all figures analyzing correlations, we colored 
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Figure 5: Above: Correlation of asymmetries S$k s versus S v >k s with the contribution of one 
mass insertion (S^ L ) 32 (left) and (Slr^^ ( r i§ht), f° r chargino (left) and gluino (right) exchanges. 
Region inside the ellipse corresponds to the allowed experimental ranges at 2a level. Below: as 
previously, but for arg^J^arg^^J (left) and arg[(<5^J 23 ]=arg[(<^ L ) 32 ] (right), with 

the contribution of two mass insertions \Sf j j i ) Sz (Sri,) (left) for gluino exchanges, and 

\ / 23 ^ * 23 

(^l_r) 23 (^ll)32 ( r ight) for both gluino and chargino exchanges. 
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Figure 6: Correlation of asymmetries S$k s versus S v ik s for chargino contribution with single 
mass insertion (5^) 32 . In the right plot the effect of a charged Higgs exchange, with mass 
rriH = 200 GeV and tan (3 = 40, has been taken into account. 



the area of the ellipse corresponding to the allowed experimental range at 2a level. 2 

In Fig. 6 the impact of a light charged Higgs in chargino exchanges is presented, when a 
charged Higgs with mass m# = 200 GeV and tan /3 — 40 has been taken into account. The 
effects of charged Higgs exchange in the case of {S^ L ) 32 mass insertion are negligible, as we 
expect from the fact that terms proportional to (S^ L ) 32 in b — * s-y and b — ► sg amplitudes 
are not enhanced by tan/?. On the other hand, in gluino exchanges with (sfj^j or 

(^ll) 23 > the most conspicuous effect of charged Higgs contribution is in populating the 
area outside the allowed experimental region. This is due to a destructive interference 
with b — > S7 amplitude, which relaxes the b — > 37 constraints. The most relevant effect 
of a charged Higgs exchange is in the scenario of chargino exchanges with (<5^ L ) 32 mass 
insertion. In this case, as can be seen from Fig. 6, a strong destructive interference 
with b — > S7 amplitude can relax the b — > 57 constraints in the right direction, allowing 
chargino predictions to fit inside the experimental region. Moreover, we have checked 
that, for tan/? = 40, charged Higgs heavier than approximately 600 GeV cannot affect 
the CP asymmetries significantly. 



4 Direct CP asymmetry in b — > versus S^^Ko 

Next we present the correlation for SUSY predictions between the direct CP asymmetry 
Acp{b —>■ sj) in b — > 57 decay and the other ones in B — > <p(r]')Ks- The CP asymmetry 

2 All ellipses here have axes of length 4cr. As a first approximation, no correlation between the expec- 
tation values of the two observablcs have been assumed. 
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Figure 7: Correlations of S^k s versus Acp(b — > S7) (left) and S^Kg versus Acp(b — > 37) (right), 
for gluino and chargino contributions with mass insertions (Slr) and (5^) 32 respectively. 

in 6 — > S7 is measured in the inclusive radiative decay of B — > X s 7 by the quantity 

- 4ot(6 - S7) = r(B-x. 7 ) + r(B-.^)- (29) 

The SM prediction for Acp(6 — > 57) is very small, less than 1% in magnitude, but known 
with high precision [7]. Indeed, inclusive decay rates of B mesons are free from large 
theoretical uncertainties since they can be reliably calculated in QCD using the OPE. 
Thus, the observation of sizeable effects in Acp{b — > 37) would be a clean signal of new 
physics. In particular, large asymmetries are expected in models with enhanced chromo- 
magnetic dipole operators, like for instance supersymmetric models [7]. 

In Fig. 7 we show our results for two mass insertions (J>lr) 23 an d {^ll)32 with both 
gluino and chargino exchanges. In this case we see that S<pK s constraints do not set any 
restriction on A C p{b — > 57), and also large and positive values of Ac pip 57) asymmetry 
can be achieved. However, by imposing the constraints on S^Kg, see plot on the right 
side of Fig. 7, the region of negative Acp(b —>■ sj) is disfavored in this scenario as well. 



5 Conclusions 

CP violation remains as one of the most interesting research topics both theoretically 
and experimentally. Especially the .B-meson decay modes seem to be ideally suited for 
searching effects of new physics, since new interesting results are mounting from the B- 
factories. For interpretation of those results, which get more accurate with more statistics, 
also reducing theoretical uncertainties in the calculations is a big challenge. 
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The 5-meson decays to <j)K, 7]'K, and to X57 provide a clean window to the physics 
beyond the SM. Here our analysis of the supersymmetric contributions to the CP asym- 
metries and the branching ratios of these processes in a model independent way has been 
reviewed. The relevant SUSY contributions in the b — > s transitions, namely chargino 
and gluino exchanges in box and penguin diagrams, have been considered by using the 
mass insertion method. 

Due to the stringent constraints from the experimental measurements of BR{b — > sj), 
the scenario with pure chargino exchanges cannot give large and negative values for CP 
asymmetry S^Ks- It is> however, seen that charged Higgs may enhance the chargino 
contributions substantially. On the other hand, it is quite possible for gluino exchanges 
to account for S^Ks an d SrfKs at the same time. 

We also discussed the correlations between the CP asymmetries of these processes and 
the direct CP asymmetry in b — > S7 decay. In this case, we found that the general trend 
of SUSY models, satisfying all the experimental constraints, favors large and positive 
contributions to b — > S7 asymmetry. 
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